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Abstract 
We take the JUNO experiment as an example to explore nuclearites in the future large liquid scintillator 
detector. Comparing to the previous calculations, the visible energy of nuclearites across the liquid scintilla- 
tor will be reestimated for the liquid scintillator based detector. Then the JUNO sensitivities to the nuclearite 
flux are presented. It is found that the JUNO projected sensitivities can be better than 7.7x 10^ "cm ?s^! sr! 
for the nuclearite mass 10? GeV < M < 107? GeV and initial velocity 107^ < Bo < 107! with a 


20 year running. Note that the JUNO will give the most stringent limits for downgoing nuclearites with 


1.6 x 10P GeV < M < 4.0 x 10? GeV and a typical galactic velocity p = 107°. 
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I. INTRODUCTION 


Strange quark matter (SQM) is a hypothetical strongly interacting matter composed of roughly 
equal numbers of u, d, s quarks and a small amount of electrons il . It is believed that SQM 
is the true ground state of quantum chromodynamics, where absolutely stable SQM objects with 
the baryon number A ranging from that of ordinary nuclei to neutron stars (A ~ 10°’) are ex- 
pected ial. The SQM has a little larger density than the saturation density of ordinary nuclear 
matter, which may be created in various situations, e.g., the hadronization process of the early 
universe [2], collision of binary compact stars |4, [5], type II supernovae driven by deconfinement 
phase transition [6], and even heavy ion collisions on Earth al. The SQM objects are consid- 
ered as the cold dark matter candidates and may be presented in the cosmic radiation reaching 
Earth. Light SQM objects (A « 107) are usually called strangelets oial 


14 lish. 


focus on heavier ones (M > 10!° GeV) known as nuclearites [14 


, while in this work we 


Based on their special 
properties, the nuclearite searches have been performed by identifying the seismic activities with 
an epilinear source on Earth and the im the ionization tracks in ancient mica and 
CR39 nuclear track detectors in the MACRO fio, bl SLIM bil. and Ohya o2 experiments, 
the bar excitations induced by the thermoacoustic effect in resonant bar detectors 3l , the Ruther- 
ford backscattering of very heavy nuclei o4 gl. the signatures of gravitational lensing caused 
by massive nuclearites dba), and the photons emitted when a nuclearite moves through water 
in the ANTARES o9) experiment and atmosphere in the future JEM-EUSO experiment, etc. 
Despite the nonobservation of nuclearites, these experiments are able to constrain the upper limits 

on the flux of cosmic nuclearites. 
The liquid scintillator (LS) as the detection medium in the past neutrino experiments has 
Ni» Now the next generation large LS detectors JUNO bd 


and LENA [37] are constructing in China and proposed in Europe, respectively. When a nucle- 


achieved great successes | 


arite passes through the LS medium, the elastic collisions between the nuclearite and ambient LS 
molecules will result in an overheating track. Many photons from the black-body radiation of this 
track can be observed by the photomultiplier tubes (PMTs). Therefore the future large LS detec- 
tors have the ability to search for nuclearites. A major advantage of the LS detectors is that the 
LS wavelength shifters can absorb the short wavelength photons and reemit the longer wavelength 
photons. This feature ensures that the LS detectors can collect more photons from the black-body 


radiation of the nuclearite track. Here we shall take the JUNO detector as an example to ex- 
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plore nuclearites. The JUNO is a 20 kton multipurpose underground LS detector and primarily 
determines the neutrino mass hierarchy by detecting reactor antineutrinos. The JUNO detector is 
deployed in a 700 m underground laboratory and consists of a central detector, a water Cherenkov 
detector and a muon tracker. The JUNO central detector holds a 20 kton LS which will be in a 
spherical container with a radius of 17.7 m Bsh. There is 1.5 m water buff region between about 
18000 20-in PMTs, 36000 3-in PMTs and the LS surface. 

In this paper, we shall explore nuclearites in the JUNO LS detector and analyze the JUNO 
detection capability. Comparing to the previous calculations, the visible energy of nuclearite per 
unit track length in the JUNO LS region will be reestimated in terms of the LS fluorescence 
quantum yields, the PMT quantum efficiencies and the JUNO detector design. Then we predict 
the JUNO sensitivities to the nuclearite flux. In Sec. II, we outline the main features of the 
nuclearite and give the maximal zenith angle below which nuclearites may pass through the Earth 
rocks and reach the JUNO detector. In Sec. III, the light yield of nuclearites traversing the JUNO 
LS will be analyzed in detail. In Sec. IV, we present the JUNO sensitivity to the nuclearite flux 


based on some conditions. Finally, some discussions and conclusions will be given in Sec. V. 


Il. THE NUCLEARITE ENERGY LOSS 


The dominant energy loss mechanism for nuclearites passing through matter is elastic or 
quasielastic collisions with the ambient atoms. As with meteorites, the nuclearite energy loss 


rate can be written as [14] 


dE 


"M 2 
dx cpp , (1) 


where £ is its velocity and p is the density of the traversed medium. The effective nuclearite cross 


section c is given by 
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where the nuclearite density is estimated to be py = 3.6x 10" e-em"? and the nuclearite radius 
Ro can be easily induced from its mass M and density py. When the nuclearite radius Ry < 1A 
(M < 8.4x10'^GeV), o is dominated by the nuclearite electron atmosphere which is never smaller 


than the typical atomic size with the radius of ~ 1 A [14 agl. It is worthwhile to stress that the 


right-hand side of Eq. should be replaced by the constant retarding force —eo for the subsonic 
velocity B < B. = Je/p with a structural energy density £ ~ 10? erg cm? i4). 

Based on Eqs. and (2), the travel length of a nuclearite depends on its mass M, velocity B 
and the medium density p. Some nuclearites may pass through the Earth rocks and arrive at the 
JUNO detector which has a 700 m rock overburden. For the JUNO detectable nuclearites, they 
will traverse different thicknesses of rocks based on their direction (zenith angle 0.). In addition, 
one should consider the change of the Earth matter density. In terms of the PREM Earth density 
profile [39], we numerically calculate the maximal zenith angle Omax below which nuclearites may 
reach the JUNO detector, namely its local velocity B; > 0 at the detector level. The corresponding 
results have been plotted in Fig. [for 10? GeV < M < 10?* GeV and five typical initial velocities 
fo at the ground level. It is clear that nuclearites from the 6, = 0° direction can reach the JUNO 
detector for M > 10P GeV and 10? < f < 107. For a typical galactic velocity By = 10^, all 
directional nuclearites can arrive at the detector when M > 2.5 x 10? GeV. Since the Earth density 


sharply changes between the core and mantle, we can see the knee points at 6, = 146.9° in Fig. 
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FIG. 1: The maximal zenith angles Omax below which the local velocity of nuclearite satisfies 8; > 0 (solid 
lines) and £4 > Bmin (dashed lines) for 5 typical initial velocities £o at the ground level. In Sec. we shall 


describe the minimal local velocity Smin- 
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Il. THE VISIBLE ENERGY OF NUCLEARITE IN JUNO LS 


When a nuclearite traverses the JUNO LS medium, the LS molecules (C;gH39) along the nucle- 
arite path will disintegrate into their constituents because of the nuclearite elastic or quasielastic 
collisions. These heated atoms will further collide with the ambient LS molecules and generate 
a hot plasma shockwave mm The evolutions of effective temperature T(t) and radius R(t) of the 


expanding thermal shockwave can be written as [14 


R(t) = V8BitRo , (3) 
T(t) = mB,Ro/(V8nt) , (4) 


where m is the LS molecule mass and n = 48 is the number of submolecular species. Here f, 
denotes the local velocity of nuclearite at the detector level. The expanding cylindrical thermal 
shockwave can emit lights through the black-body radiation. The corresponding power spectrum 
is given by 


dp ho? 1 
doda — 4c? eholkr —] ’ 


(5) 
where w is the angular frequency and a denotes the area of shockwave. Then we can deduce the 
emitted photon numbers per unit track length dN, /dx from the expanding cylindrical shock: 


dN. 1 
y f di Í diano- A 


dx dwda w 


8: ig l 
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Tm 


where we have used the natural system of units with A = c = k = 1. Similarly, the total emitted 
energy dE,/dx can be directly obtained through the replacement w? —> œ? in Eq. (6). 

These emitted photons from the black-body radiation cannot be entirely detected by the JUNO 
PMTs since they will suffer the absorption, reemission and Rayleigh scattering processes in the 
JUNO LS [36]. On the other hand, the PMT quantum efficiency is related to the photon wave- 
length. Therefore one cannot simply use the total emitted energy dE, /dx or the total photon num- 
bers to describe the visible energy in the LS detector. It is convenient for us to calculate the visible 
energy if the photon electron (pe) efficiency per photon e(4) is available for the JUNO detector. 
Based on the LAB, PPO and bis-MSB fluorescence quantum yields [40], we adopt a combined 
PMT quantum efficiency curve shape from the Hamamatsu PMT data (400-800 nm) and a 
fixed 27% efficiency (250-400 nm) to calculate e(4) for the wavelength range 250 nm < A < 800 


5 


nm as shown in Fig. We have assumed an averaged 60% survival probability of reemitted 
photons from the detector center to the PMT surface and the 75% PMT photocathode coverage. It 
is found that the modeled e(A) approaches to zero for A > 640 nm. In the absence of the related 
experimental data, we do not include the contribution of the 2 < 250 nm case. Note that it will not 


significantly affect our final results. 
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FIG. 2: The modeled photon electron efficiency per photon as a function of the wavelength 4. 


With the help of Eq. (6) and the photon electron efficiency e in Fig. 2] we can deduce the 


visible energy of nuclearite per unit track length in the JUNO LS: 
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where fmin takes the larger one of t) = Ro/( V861) and t, = (1/Ro)*to with the mean free 
path / ~ 2.7A. In Eq. (7), we have simply used that the 1 MeV gamma in the detector center 
may averagely produce 1200 photon electrons m Based on Eq. (7), one may numerically 
calculate dE,;;/dx as shown in the left panel of Fig. It is clear that dE,;;/dx does not vary for 
M < 84 x 10^ GeV. This is because we have adopted a constant value for Ro in Eqs. (3) and (4), 
i.e., the radius of nuclearite electron atmosphere 1 A. In the right panel of Fig.|3| we plot the ratio 
of dE, /dx and dE/dx as a function of the local nuclearite velocity B,. It is found that this ratio is 
independent of the nuclearite mass M when M < 8.4 x 10^ GeV and M > 1.7 x 10!6 GeV. For 
the M > 1.7 x 10/6 GeV case, namely fmin = fo > tı, one may easily find dE, /dx œ R? through 


the variable substitution t = t’Ro. 
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FIG. 3: Left panel: the visible energy of nuclearite per unit track length dEy;; /dx as a function of nuclearite 


mass for 5 typical 61; Right panel: the ratio of dEyis/dx and dE/dx as a function of Bj. 
IV. THE EXPECTED JUNO SENSITIVITIES 


The light signals from nuclearites can be recorded when they satisfy the JUNO trigger condi- 
tions. Here we assume that the JUNO trigger threshold is 0.5 MeV within a 300 ns window for 


the following analyses. Then one may obtain 
dE, 
f x 300ns x ze 20.5 MeV. (8) 
x 


With the help of Eqs. (7) and (8), we calculate the minimal local velocity Bmin as shown in Fig. 
It is clear that the local velocity £; must be larger than 8.7 x 10~° for 10? GeV < M x 10? GeV. 
For a fixed initial velocity o at the ground level, the maximal zenith angle Omax can be deduced 
from Eq. and the requirement 8; > min. In Fig. [1] we have plotted the corresponding Omax 
with dashed lines. It is found that the JUNO may detect all downgoing nuclearites (zenith angle 
6. < 90°) with 5.0 x 10? GeV < M x 10? GeV and 10? < f < 107!. For the Bọ = 107 case, 
the downgoing nuclearites with M > 6.3 x 10” GeV can satisfy the 61 > Bmin condition. 


The expected nuclearite numbers in JUNO can be written as 


Ns = 2n (1 — cos Omax) Ø Trun 7 Rog (9) 


where Tmn is the JUNO running time and ¢ is the isotropic nuclearite flux in unit of cm ?s !sr !. 


Here we require that the nuclearite track length in LS region should be larger than 5 m and derive 


the effective JUNO radius Reg = 4(07.7my — (5m/2)? = 17.52 m. Then the 90% confidence level 
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FIG. 4: The minimal local velocity Smin above which nuclearites can satisfy the assumed JUNO trigger. 


(C.L.) upper limit Noo to the expected Ns can be derived through the following formula lla. 


f°, Lass INs)dNs 
90% = =E (10) 
Jio LNovslNs Ns 
with the Poisson-based likelihood function 
Ns + Nac) 
L(Nops|Ns) = 2S + Na) caeno (11) 


Nobs! 


To estimate the JUNO sensitivities to the nuclearite flux ¢, we assume the background number 
Nac = 0 and take the observed event number Nos, = Nag = 0 for a 20 yr running. With the help 
of Eqs. (9) and (10), we plot the 90% C.L. flux upper limits (solid lines) for five typical initial 
velocities By as shown in the left panel of Fig. [5| It is clear that the JUNO sensitivities are better 
than 7.7 x 107" cm^s"!sr^! for 10^ GeV x M x 10% GeV and 107* x Bp x 107!. The JUNO is 
only sensitive to a narrow parameter space of the By = 10^? case because of B, < Bmin. In addition, 
the most optimistic limit (black dotted line) has been also plotted for the By = 10? case where we 
only require 6, > 0 and take Reg = 17.7 m. 

To compare with the MACRO bd. ANTARES bd. SLIM and Ohya experimental 
results, we calculate the JUNO upper limit on the downgoing nuclearites in the By = 10? case. Our 
numerical results are presented in the right panel of Fig. [5| It is clear that the JUNO sensitivity 
is far better than the MACRO, SLIM and Ohya limits. Note that the JUNO will give the most 


stringent limit in the range of 1.6 x 10? GeV < M x 4.0 x 10? GeV. In the most optimistic case 
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FIG. 5: The JUNO 90% C.L. flux upper limits to the all direction (left) and downgoing (right) nuclearites 
with the local velocity £1 > Smin and the effective LS radius Reg = 17.52 m for a 20 yr running. The black 
and red dotted lines describe the JUNO most optimistic limits with 8; > 0 and Reg = 17.7 m in the initial 


velocity Bo = 10^? case. 


(red dotted line), the above range can be extended to 5.0 x 10? GeV «x M x 4.0 x 10P GeV. 
Here we have also plotted the galactic dark matter (DM) upper limit $4; = opwo/(21M) 
where nuclearites are assumed to contribute all of the local DM density ppm = 0.39 GeV cm? 
ld. For M > 3.1 x 10” GeV, the galactic DM limit is dominant. In the future, the JEM-EUSO 


experiment will give a more stringent limit ¢ < 10? cm ?s !sr ! for M > 5x 10? GeV 30). 


V. DISCUSSIONS AND CONCLUSIONS 


As mentioned in Sec. the photon electron efficiency per photon e(A) is not considered for 
the A « 250 nm range because of the absence of the related experimental data. If these data are 
available in future, we shall derive the larger d E,;,/dx than those in the left panel of Fig. B| Then 
the smaller Smin will be expected. The predicted sensitivities (solid lines from 8, > min) in Fig. 
will approach the most optimistic limits (the corresponding dotted lines from 6; > 0). It is clear 
that our results do not change significantly for M > 1.0 x 10^ GeV and By = 107°. Note that 
the nuclearite mass cannot be correctly reconstructed from the incomplete e(4) when a nuclearite 
is really detected by the JUNO LS detector. In addition, the JUNO can only give the mass lower 
bound for very large dE,;,/dx because of the PMT saturation. 

In conclusion, we have investigated nuclearites in the JUNO LS detector. Comparing to the 


previous calculations, the visible energy of nuclearite in the LS has been estimated in detail. Then 
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we give the JUNO detectable range of the zenith angle for the nuclearite mass 10? GeV < M < 
10°% GeV and five typical initial velocities By at the ground level. Finally, we present the JUNO 
sensitivities to the nuclearite flux for a 20 yr running. It is found that the JUNO sensitivities to 
all directional nuclearites are better than 7.7 x 10 "cm ?s !sr ! for 10'° GeV < M x 10% GeV 
and 10^ < By x 10^!. For the downgoing nuclearites, the expected sensitivities are much better 
than those from the MACRO, SLIM and Ohya experiments in the case of 8) = 1073. Note that the 
JUNO will give the most stringent limits for 1.6 x 10? GeV < M x 4.0 x 10? GeV. 
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